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Abstract: 

Models with a dynamic cosmological term A {t) are becoming popular as 
they solve the cosmological constant problem in a natural way. Instead of 
considering any ad-hoc assumption for the variation of A, we consider a 
particular symmetry, the contracted Ricci-coUineation along the fluid flow, 
in Einstein's theory. We show that apart from having interesting properties, 
this symmetry does demand A to be a function of the scale factor of the 
Robertson- Walker metric. In order to test the consistency of the resulting 
model with observations, we study the magnitude-redshift relation for the 
type la supernovae data from Perlmutter et al. The data fit the model very 
well and require a positive non-zero A and a negative deceleration parameter. 
The best-fitting fiat model is obtained as Qq ~ 0.5 with go ~ —0.2. 
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1. INTRODUCTION 

The recent measurements of the CMB anisotropy [1] and the observations of 
type la supernovae made, independently, by Perlmutter et al [2] and Riess 
et al [3] consistently demand a significant and positive cosmological constant 
A. These observations suggest Friedmann models with negative pressure- 
matter in which the expansion of the universe is accelerating. Observations 
of gravitational lensing [4] also indicate the presence of a non-zero A. 

On the other hand, a dynamical A(t) has been considered in numerous 
papers in order to explain its observed small value which is about 120 orders 

iJUCAA, Post Bag 4, Ganeshkhind, Pune 411 007, India 
E-mail: vishwa@iucaa.ernet.in 



1 



of magnitude below the value for the vacuum energy density predicted by 
quantum field theory- the so-called cosmological constant problem [5]. This 
phenomenological solution is based on the argument that A relaxed to its 
present estimate due to the expansion of the universe. (It is customary to 
associate a positive cosmological constant A with a vacuum energy density 
Pv = A/SttG.) 

As the dynamics of the variable A-models depends sensitively on the cho- 
sen dynamic law for the variation of A and, in general, becomes altogether 
different from the dynamics of the corresponding constant A-models, there is 
no reason to believe that the observations of distant objects would also agree 
with the variable A-models, given that they agree with the corresponding 
constant A-ones, especially for the same estimates of the parameters. In this 
view, it would be worth while to test the consistency of the above-mentioned 
observations with the variable A-models which solve the cosmological con- 
stant problem in a natural way. It may be mentioned that the quintessence 
or the exotic 'x-fluid' models (— 1 < Px/ Px < 0) are also capable to explain, in 
a natural way, the smallness of the present vacuum energy density, since the 
associated vacuum energy density dynamically evolves towards zero with the 
expansion of the universe. However, as Garnavich et al [6] have shown, if one 
attempts to constrain the equation of state of the dark energy component 
(that may have contributed to accelerating the cosmic expansion), with the 
la supernova data, one finds that the data is consistent with either a cosmo- 
logical constant, or a scalar field that has, on average, an equation-of-state 
parameter similar to the one for the cosmological constant, i.e., -1. However, 
we note that Px/Px = ~1 corresponds to the cosmological A which, in view 
of the conservation of the energy momentum tensor of the exotic component, 
becomes a true constant and, hence, cannot solve the cosmological constant 
problem. However, as we shall see later, the usual kinematical A can always 
be, in general, a function of time as the conserved quantity in this case is 
matter plus vacuum, and not the vacuum only [7]. 

We note that a number of models with A as a function of time have 
been presented in recent years [8, 10-13]. Different phenomenological laws, 
for the decay of A, have been proposed in these models which are either 
from dimensional arguments or from ad-hoc assumptions. Though the pre- 
cise mechanism of vacuum-decay, which should come from the fundamental 
interactions, is not yet known, it would, however, be worth while to look for 
some symmetry principles which actually demand the variation of A. More- 
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over, it is always reasonable to consider symmetry properties of spacetime 
rather than considering ad-hoc assumptions for the variation of A. 

In this paper, we consider, on the level of classical general relativity, a 
particular Ricci-symmetry which is the contracted Ricci-coUineation along 
the fluid flow vector and show that, apart from having interesting properties, 
this symmetry does demand A to be a function of time (and space, in general). 
This is done in section 2. In section 3, we derive the magnitude (m) - 
redshift (z) relation in the resulting model to test its consistency with the la 
supernovae data from Perlmutter et al [2] . The model fits the data very well. 
The numerical results are discussed in section 4 followed by a concluding 
section. 

Earlier some cosmologicai models with this symmetry were discussed by 
considering particular sets of initial and boundary conditions and by assum- 
ing that A = A{t) in some of them [8-10]. However, as wc shall see in section 
2, the incorporation of this particular symmetry in the Einstein theory does 
demand A to be a variable and there is no need to assume A = A(t) a priori. 
For the ready reference and completeness, we reproduce the model in the 
following. 



2. SYMMETRY CONSIDERATIONS AND THE RESULTING 
MODEL 

We begin our discussion by considering the contracted Ricci-coUineation. 
The motivation for considering this particular symmetry will be discussed 
later. We note that a spacetime is said to admit a Ricci-coUineation along a 
field vector 77' if [14] 

Cr]Rij — 0, (1) 

where denotes the Lie-derivative along 77^. Further, a spacetime is said to 
admit a family of contracted Ricci-coUineation if 

g'^Cr,Rij = 0, (2) 
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which leads to the conservation law generator^ 

if the Einstein field equations 

R^i - iRg^^ = -8nG (t^^ - ^g'^^ (4) 

are satisfied. We consider units with c = 1. Recalling that the energy density 
associated with vacuum is Pv = A/SttG with its pressure = — Pv implying 
T^^ = —{A/87rG)g^\ the quantity appearing in the parentheses on the right 
hand side of equation (4) can be written as 

Tg = + Tl^ = (ptot + Vtot)u'u^ + Vtotg''. (5) 

which represents the energy momentum tensor of the total matter, i.e., or- 
dinary matter plus vacuum. Here ptot = P + Pv, Ptot = P + Pv = P — Pv and 

is the normalized velocity 4-vector. The Bianchi identities then confirm 
the conservation of T^t and not the conservation of T*-' and T^^ separately 
unless A = constant. 

If we consider rj^ oc u\ equation (3) reduces to 

{(Ptot + 3PtotW};j = 0, (6) 

which may be interpreted as the conservation of generalized momentum den- 
sity. This is an important result in its own right , but it implies even more. 
To understand the full meaning of this conservation law and compare it with 
the existing results, let us consider the Robertson- Walker metric 

ds^ = -dt" + S\t) + r\de'' + sin2 edct)'')] , (7) 

[ 1 — kr^ J 

characterized by the scale factor S{t) and the curvature index A;g{— 1,0,1} 
of the spatial hypersurfaces t = constant. Now the Einstein field equations 
(4) yield two independent equations: 

^Due to a typographical error in [8], the same definition of the conservation law gen- 
erator for a contracted Ricci-coUineation was also assigned to the Ricci-collineation. But 
this did not affect the results of the paper as the Ricci-collineation was never used in the 
paper. 
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the Raychaudhuri equation: 



S 
S 



AttG 
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(Ptot + 3ptot) 



(8) 



and the Priedmann equation: 



k SttG 



(9) 




By using (7), the conservation law (6) reduces to 



(Ptot + 3ptot)'S'^ = constant A (say), 



(10) 



which is the central equation of our analysis and may be interpreted as the 
conservation of total active gravitational mass, taken matter and vacuum 
together, of a comoving sphere of radius S. Note that (p + 3p) is defined 
as the active gravitational mass density of the universe [15]. Obviously the 
conservation law (10) reduces to the pressure-less phase of the standard big 
bang (FLRW) model for p = A = 0. 

To understand the presence of pressure and vacuum terms in equation 
(10), we consider the Raychaudhuri equation (8) which, being the analogue 
of Newtonian gravitation, suggests that the gravitational attraction, experi- 
enced by a unit test mass, is exerted in fact not only by p as in the Newtonian 
theory but rather by the active gravitational mass density (p + 3p), which 
exhibits the relativistic effects. Associated with the attractive force GA/S'^ 
{A= ^{p + 3^)5"^ being the active gravitational mass from the gravitating 
matter) in equation (8), there is a repulsive force — A5'/3 due to a positive A 
(or an additional attraction with a negative A). It is this repulsive force which 
drives inflationary expansion in the early vacuum-dominated universe. It is 
obvious that the constant A (the total active gravitational mass) might take 
different values in different phases of evolution depending upon the relative 
dominance of these two terms of opposite character in equation (8). This im- 
plies that the total active gravitational mass might be conserved phase- wise 
only and not throughout the evolution. (For details, see [8]). 

One may also note that the left hand side of equation (8) is the Gaussian 
curvature of the two-dimensional surface specified by varying r and t, keeping 
9 and constant, in equation (7) and may be considered as the curvature 
of the homogeneous and isotropic spacetime. Thus equation (8) implies that 



5 



the curvature of spacetime is governed by the total active gravitational mass 
density of the universe. In this view, a naive assumption would be that the 
cause of curvature of spacetime, i.e., the total active gravitational mass be 
conserved, justifying our symmetry consideration leading to the conservation 
law (10). Equation (10) also implies, via equation (8), that the curvature of 
spacetime evolves as 5*^'^, which quickly transforms the spacetime from a 
state of large curvature to a state of flatness. 

It would be worth while to mention that if one wishes to describe the 
radiation-dominated universe in the Newtonian framework (extending the 
work of McCrea and Milne), one needs to use (p + 3p) instead of p only in 
Poisson's equation to get the right answer. This is consistent with our using 
(p + 3p) for gravitational mass density. 

Equations (8) and (10), taken together, yield 



which integrates to 



S^^^^-rB, (12) 

where 5 is a constant of integration. This supplies the dynamics of the scale 
factor. Equations (9), (10) and (12) may be used to obtain 

A 3(3 + k) 
Ptot - ^3 + SnGS^ ' ^^"^^ 

(B + k) 

" "^^i^' ^^^^ 
which give the equation of state of the perfect fluid constituting the total 
matter of the universe as 



Pt'ot = ^(Ptot + 3ptot)', (15) 

where K — —{B + k)^/ {S7rG)^A'^. This is a physically reasonable equation of 
state since dptot/dptot = 2ptot/3(ptot +Ptot) indicating that dptot/dptot < 1/3 
for Ptot > Ptot (provided (ptot +Ptot) > 0). It may be noted that the equation 
of state (15) breaks down for ptot = (in the same way as does the usual 
barotropic equation of state for a perfect fluid for p = 0), in which case 
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equations (13) and (14) may be treated as the parametric equations of state. 
Consequences of the resulting models for the case A = have been discussed 
elsewhere by Abdussattar and Vishwakarma [9] where the models obviously 
get constrained hj A > and {B + k) < and by the additional constraint 
S < 8ttGA/3 \B + k\ when {B + k) < 0. 

When A 7^ 0, the equation of state (15) docs not supply information on 
the ordinary matter source. If the ordinary matter is baryonic with its usual 
barotropic equation of state 

p — wp, < -u; < 1, (16) 

this simply implies that A cannot just remain constant. The reason is obvi- 
ous. We now have 4 independent equations (12)-(14) and (16) in 3 unknowns 
S, p and p. This over-determinacy can be compensated by allowing at least 
one of the remaining parameters to vary. The only such parameter of interest 
is A if we keep G as a constant (the cases with variable G have been discussed 
elsewhere [10]). It is obvious that we would have reached the same conclu- 
sion, had we considered any other assumption in place of (2). However, this 
symmetry, as we have seen, has its own significance. Equations (13), (14) 
and (16) thus yield 



A 



SttG 



wA ^ (l + 3w)(5 + fc) 



1 



53 MGS^ 
'A (B + k) 



'l+w) 



(17) 



(18) 



53 AttGS^ 

It may be mentioned that A varying as S*^^, which the present model 
has in the pressure-less phase of evolution, has also been considered by sev- 
eral authors to explain the present small value of A [11, 12]. The ansatz is 
primarily due to Chen and Wu [11] who postulated it through dimensional 
arguments made in the spirit of quantum cosmology. Recently Jafarizadeh 
et al [16] have calculated the tunneling rate with a cosmological constant de- 
caying as S^"^ and concluded that the most probable cosmological term with 
the highest tunneling rate occurred at m = 2. However, the present model 
differs from the above-mentioned models in the sense that contrary to the 
A ~ S~'^ throughout the evolution (as has been assumed in these models), A 
varies differently in the different phases of evolution in the present model as 
is clear from equation (17). 
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While comparing the model with SN la data, we shall be interested in 
the effects which occurred at redshift < 1. We, therefore, neglect radiation 
and consider w — in the following. (The parameters in the early phase of 
evolution can be calculated by following [8]). With this in view, we recast 
equations (8) and (9) in the following forms to give the relative contributions 
of the different cosmological parameters at the present epoch: 

2 [go + f^Ao] = ^^0, (19) 

1 + QfeO = ^^0 + ^^AO, (20) 

where Q = SnGp/SH'^, Qk = k/S'^H'^ and Qa = A/3if^ are, respectively, the 
dimensionless forms of the density, the curvature and the cosmological con- 
stant parameters, and the subscript characterizes the value of the quantity 
at the present epoch. Equations (10) and (12) may be used to give the values 
of the constants A and B as 

0C3O-2 

A=(Qo-2Qao)^^|^, (21) 
B = (2Qao - Qo + l)HlSl. (22) 



3. MAGNITUDE-REDSHIFT RELATION IN THE MODEL 

The cosmic distance measures, like the luminosity distance and the an- 
gular size distance^ depend sensitively on the spatial curvature and the ex- 
pansion dynamics of the models and consequently on the present densities of 
the various energy components and their equations of state. For this reason, 
the magnitude-redshift relation for distant standard candles and the angular 
size-redshift relation for distant standard measuring rods have been proposed 
as potential tests for cosmological models and play crucial role in determining 
cosmological parameters. 

Let us consider that the observer at r = and t — tQ receives light emitted 
aX t — ti from a source of absolute luminosity L located at a radial distance 
ri. The cosmological redshift z of the source is related with ti and to by 
\ + z = S(to)/ S(ti). If the (apparent) luminosity of the source measured by 
the observer is /, the luminosity distance d]^ of the source, defined by 
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is then given by 



(23) 



dj^ = {l + z)Son. (24) 

For historical reasons, the observed and absolute luminosities / and L are 
defined, respectively, in terms of the K-corrected observed and absolute 
magnitudes m and M as / = 10^^™/^ x 2.52 x 10^^ erg cm~^ and 
L = 10-2A^/5 X 3.02 X 10^^ erg s"^ [17]. When written in terms of m and M, 
equation (23) yields 

m{z; M, Oo, ^Ao) = M + 5\og,Q{Vj^{z; Qq, I^aq)}, (25) 

where M = M — 51ogio-f^o + 25 and Vi^^z; fio, I^ao) = Hq di^{z; Qq, ^l^o, Hq) 
is the dimensionless luminosity distance. Here dj^ is measured in Mpc. By 
using equation (7), the coordinate distance ri, appearing in equation (24), 
yields 

V'(ri) = r ^ (26) 

JSo/il+z) SS 

with 

■0(ri) = sin"''^ ri, k — 1 
= ri, k = 

= sinh"Vi, k = -l. (27) 
By the use of (12), (21) and (22), equation (26) yields 

= ^ r[{'2^Ao-^o + m+^T-i'^^Ao-^o){l+z')T'^''dz'. (28) 

Equations (24), (27) and (28) can also be combined into a single compact 
equation as 
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e (v^ _^'[(2fiAo - 1^0 + 1)(1 + z')' - (2fiAo - Qo){l + z'f]-'/'dz'^ , 

(29) 

where 

^(x) = sin(a;) with /C = Qko when Jlfeo > 0, 
^{x) = sinh(x) with /C = —Qko when < 0, 

^(x) = X with /C = 1 when Q^o = 0. 

Thus for given Ai, flo and f^Ao, equations (25) and (29) give the predicted 
value of m{z) at a given z. By using the K-corrected effective magnitudes 
mf^, which have also been corrected for the light-curve width-luminosity 
relation and galactic extinction, and using the same standard errors az,i and 
(7 cflF of the ith supernova with redshift Zi as used by Perlmutter et al, we 

i 

compute according to 



The best fit parameters are obtained by minimising this equation. We note 
that equation (29) is sensitive to Qq and Qao for distant sources only. For 
the nearby sources (in the low-redshift hmit), equations (25) and (29) reduce 
to 

m{z) = M + 5 logioz, (31) 

which can be used to measure A4 by using low-redshift supernovae- measurements 
that are far enough into the Hubble flow so that their pecuhar velocities do 
not contribute significantly to their redshifts. 



4. NUMERICAL RESULTS 

We consider the low-redshift data on m^^ and z from the Calan-Tololo 
sample of 16 supernovae (excluding 2 outliers from the full sample of 18 super- 
novae) to estimate M. from equations [30] and (31). This gives M. — 24.03 
(in units with c = 1). This value is used in equation (25) to estimate flo 
and Qao from the high-redshift measurements. For this purpose we consider 
the data set on rrf^ and z of 38 supernovae from the Supernova Cosmology 
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Project (excluding 2 outliers and 2 likely reddened ones from the full sam- 
ple of 42 supernovae). We perform a two-parameter fit of this data set by 
following the fitting procedure of Perlmutter et alQ considered in their Fit M 
(see Fig. 5 (f) in their paper [2]). 

We find that the data shows a very good fit to the model giving the 
best-fitting flat model (fi + = 1) as = 0.54 with = 45.34, which 
shows that the fit is almost as good as to the constant A-fiat FLRW model: 
fio = 0.40 with = 44.92. 

The minimisation process gives the global best-fitting solution (calculated 
by giving free rein to Qq and f^Ao) as Qq = 1-76 and f^Ao = 1-34 with = 
44.78 which shows a rather high value of Qq and does not seem realistic in 
view of the small observed value of Qo- One can see that the data predicts an 
accelerating expansion (go < 0) of the universe as in the constant A-FLRW 
model. 

The fit of the fiat model to the actual data points has been shown in 
Figure 1, where we have compared it with the simplest A = model, i.e., 
the Einstein-de Sitter model {Q = 1), which is ruled out by the data (x^ = 
115.81). 

It may be mentioned that elsewhere [7], we have used this data in the 
model, with a different fitting procedure, by fitting the low- and high-redshift 
measurements simultaneously to equation (30). The best-fitting solutions so 
obtained are in good agreement with those obtained here. 

It would also be worth while to mention that the present model is con- 
sistent not only with the supernovae-data but it also fits very well the data 
on the angular size and redshift of the ultracompact radio sources complied 
by Jackson and Dodgson [19] as well as the updated and modified data on 
compact radio sources from Gurvits et al [7]. 

^ It may be mentioned that Perlmutter et al [2] have also fitted the data for only 3 
parameters , fio and flAo and not for the 4 parameters a, M, flo and Qao as mentioned 
in their paper (from a personal discussion with Professor R. S. Ellis, one of the authors of 
the paper). A self-consistent 4-parameter fit has been done by Efstathiou et al[18]. 
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Figure 1. Hubble diagram for 38 high-redshift and 16 low-redshift su- 
pernovae: The solid curve represents the best-fitting fiat model {Qq = 
0.54, fl\o = 0.46). For comparison, the canonical Einstein-de Sitter model 
{Q.k = 0, r^A = 0) has also been plotted (dashed curve). 



4. CONCLUSION 

In order to solve the cosmological constant problem, which has been made 
even more acute by the consequences of the current observations of CMB and 
type la supernovae, several authors have invoked a variable cosmological term 
A{t). Instead of considering any ad-hoc assumption for the variation of A, as 
has been mainly done by the authors, we have considered a particular sym- 
metry of the homogeneous, isotropic spacetime which is the contracted Ricci- 
collineation along the fiuid fiow. It have been shown that the incorporation 
of this additional symmetry into the RW model filled with baryonic matter 
renders the cosmological term A a decaying function of the scale factor. This 
helps in solving the cosmological constant problem. This new symmetry in 
the model leads to the conservation of the total active gravitational mass of 
the universe. 



12 



The resulting model fits the SN la data from Perlmutter et al very well 
and requires an accelerated expansion of the universe with a non-zero positive 
cosmological term. The goodness of fit of the data to the model is almost 
the same as in the case of the constant A-FLRW model. The estimates of 
the parameters for the best-fitting fiat model are obtained as Qq = 0.54 
and Qao = 0.46. However, the global best-fitting solution, Qq = 1.76 with 
= 1.34, does not seem realistic (as is the case with the constant A-FLRW 
model; see reference [2]) in view of the small observed value of Qq. It is also 
noted that the estimates of the density parameter Qq for the present variable 
A-model are found a bit higher than those for the constant A-FLRW model. 
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